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We measure two higher-order power spectra involving weighted cubic and squared temperature 
anisotropy maps from WMAP 5-year data to study the trispectrum generated by primordial non- 
Gaussianity. Using these measurements combined with Gaussian and noise simulations, we constrain 
the cubic order non-Gaussianity parameters tnl, and 3nl- With V-(-W-band data out to /max = 600, 
we find —7.4 < (jinl/IO^ < 8.2 and —0.6 < tnl/IO* < 3.3 improving the previous COBE-based limit 
on Tnl < 10* nearly four orders of magnitude with WMAP. We find that the ratio of trispectrum 
to bispectrum amplitude as captured by the ratio of tni/(6/nl/5)^ ranges from -3 to 21 at the 95% 
confidence level. 



Introduction. — Cosmic Inflation has deservedly be- 
come a cornerstone of modern cosmology [IH3] . Inflation 
solves the flatness, horizon and the monopole problems 
of the standard Big-Bang cosmology. Furthermore, in- 
flation is the prevailing paradigm related to the origin 
of density perturbations that gave rise to the large-scale 
structure we see in the universe. It posits that a nearly 
exponential expansion stretched space in the first mo- 
ments of the early universe and promoted microscopic 
quantum fluctuations to perturbations on cosmological 
scales today [H |3] . 

In the simple scenario, inflation is driven by a single 
scalar field whose potential energy dominates it's kinetic 
energy. This "slow-roll" situation leads to an exponen- 
tial expansion of the cosmic spacetime via the Einstein 
equations coupled to a scalar field. In order to maintain 
a slow roll, the scalar field must have minimal self inter- 
actions. Such a non-interacting field has the statistical 
feature that its fluctuations are Gaussian. A Gaussian 
fleld contains no correlations between points arising from 
interactions. 

In addition to the standard inflationary scenario, sev- 
eral possible mechanisms of inflation have emerged. 
These models usually involve multiple fields and or 
other exotic objects such as branes (motivated by string 
theory) that have non-trivial interactions. These in- 
teractions produce a departure from Gaussianity in 
a model-dependent manner [BHS]. Constraining non- 
Gaussianity therefore is important to distinguish between 
the plethora of inflationary models [lO] . 

The first order non-Gaussian parameter, /nl; has been 
measured with increasing success in the bispectrum, or 
three-point correlation function of temperature anisotr- 
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pies of the cosmic microwave background (CMB). Such 
studies have found /nl to be consistent with zero pTHTl] . 

In addition to /nl, with the trispectrum or four point 
correlation function of CMB anisotrpies, we can mea- 
sure second order non-Gaussian parameters tnl and ^nl- 
Constraints on tnl have never been directly preformed 
with data. The very weak constraint often quoted in the 
literature with tnl < 10* is an estimate based on a null 
detection of the COBE data [13 US]. Furthermore, 
the ratio between tnl and (6/nl/5)^, involving the am- 
plitudes of trispectrum and bispectrum, is often is often a 
constraint on inflationary models [51 [S]. In fact, for many 
inflationary models this ratio is related to the tcnsor-to- 
scalar ratio, which is now probed by CMB polarization 
data, providing additional consistency relations to under- 
stand the underlying physics of inflation [9' . Constraints 
on gNL have only recently been preformed, but not with 
the trispectrum of the CMB [TT] . 

In this analysis we use the trispectrum, or the four- 
point correlation function of temperature anisoptropies, 
to measure primordial non-Gaussianity at second order 
using WMAP 5-year data [T8] . 

Theory. — To parameterize the non-Gaussianity of a 
nearly Gaussian field, such as the primordial curvature 
perturbations <I'(x), we can expand it perturbatively [15] 
to second order as: 

$(x) = 0l(x) + /nl [</'i(x) - (0i(x))] +gNL0i(x) (1) 

where 0l(x) is the purely Gaussian part with /nl and 
f/NL parametrizing the first and second order deviations 
from Gaussianity. Fortunately, information about the 
curvature perturbations are contained within the CMB 
through the spherical harmonic coefficients of the tem- 



2 



perature anisotropies: 



manner Ci is the angular power spectrum of temperature 
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where <&(k) are the primordial curvature perturbations, 
gTi is the radiation transfer function, 9 is the field of 
temperature fluctuations in the CMB and Y^'s are the 
spherical harmonics. 

If the curvature perturbations are purely Gaussian, 
all the statistical information we can say about them 
is contained in the two point correlation function 
(<i>(xi)$(x2)). The information contained in the two 
point function is usually extracted in spherical harmonic 
space, leading to the power spectrum C;, defined by: 
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However, if the curvature perturbations are slightly non- 
Gaussian, this two point function is no longer suffi- 
cient to articulate all the information contained in the 
field. With non-Gaussianity, extra information can be 
extracted from the three, four and higher n-point corre- 
lation functions. 

In this paper we look to the trispectrum, or four 
point correlation function, to probe non-Gaussianity. 
To extract information from the four point function 
(<I>(xi)$(x2)<I'(x3)<I>(x4)), we again work in spherical 
harmonic space and compute: 

(5) 

(^/imi ^/2^'^2 ^^3^^3 ^^4^4 )g ^ ^^2^2 ^^3^i3 ^^4'^4 ) C ' 

where we see that the four point function breaks up into a 
piece representing a Gaussian contribution and a second 
piece representing the non-Gaussian contribution. The 
non-Gaussian piece is commonly referred to as the con- 
nected piece and is generated by interactions between 
points. 

As with the two point function, from the four point 
function we can create various power spectra. Two in- 
dependent spectra that can be produced are IC'f'^^ and 
/Cp'^^ (See Fig. [l| that are related to various compres- 
sions of ajmS as follows: 
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Here ^ should be read as is related to as the full expres- 
sion is very complex. (We refer the reader to Refs. [20] for 
full details of the derivation of eq. [6] including the pre- 
factors.). The superscript in the a|^^refers to how many 
temperature maps must be combined to form the aim- 
For example, only one is used in eq. 3]|4 In the same 



and /Cp'^'' and 1C\'^''^' are the angular power spectrum 
of squared-temperature correlated against squared tem- 
perature and cubic temperature correlated against tem- 
perature, respectively. To optimize for primordial non- 
Gaussianity detection, however, the temperature maps 
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used for /Cp'^' and 1C\^'^' are weighted appropriately in 

statistic for the 
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the same manner we optimized 
bispectrum measurement |14j. 

The estimators /Cp'^'' and /Cp'^'' contain contributions 
from the Gaussian and connected piece of equation [5] 
The connected piece further breaks up into two indepen- 
dent pieces parameterized by two parameters tnl and 
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Here, A^^'""^ and Bp'''^ are the theoretical expectations 
for the connected part of /Cp'^'' and /Cp'^-* with the as- 
sumption that Tnl and (^nl = 1 respectively. 

It is clear from equation [T] that ^nl is a non-Gaussian 
parameter of second order. Additionally, in a derivation 
of the trispectrum one finds that tnl is related to 
in a model dependent way, and is therefore also a second 
order parameter. 

Since the exact relation between tnl from the trispec- 
trum and /nl from the bispectrum is model dependent, 
this relationship gives another test to distinguish between 
models for primordial perturbations [BHS]. It is common 
to compare tnl with (6/nl/5)^ and therefore we con- 
strain the relation Anl = tnl/(6/nl/5)^ where Anl is 
the amplitude of the trispectrum to the squared bispec- 
trum. For many models Anl can be used to further con- 
strain the tensor-to-scalar ratio r |0. 

Analysis and Results. — Our recipe for analysis is 

1. We calculate Af^'^\ 6p'^\ ^p'^' and Sp'^^ in 
Eq. [?]■ [Sjfor Tnl and ^nl — 1- (See [5D] for more 
details on how to do this and other steps.) 

2. We extract /Cp'^^ and /Cp'^^ directly from WMAP 
5-year data. 

fSl) (22) 

3. We preform the extraction of /C; ' and /Cj ' ' from 
250 Gaussian maps, allowing us to determine error 
bars and the Gaussian piece of each estimator. 

4. We subtract off the Gaussian contribution to these 
estimators to ensure we are fitting to the non- 
Gaussian contribution. 

5. We fit the two unknowns tnl and (Jnl from data 
using the two equations simultaneously. The am- 
plitudes the theoretical curves must be scaled by 
gives the values for tnl and c/nl 
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FIG. 1: The top plot shows the /Cj ' and /C; ' estimators, 
shown in green and blue respectively, taken from data for the 
W band. The same estimators for the V band are shown on 
the bottom. Additionally on the top the theoretical contri- 
butions for /Cf'^ and ACf'^ proportional to tnl are shown with 
the bottom showing those proportional to (?nl. The Gaussian 
contributions were not removed from these plots. 



6. We constrain Anl by comparing tnl from the 
trispectrum with (6/nl/5)^ coming from the bis- 
pectrum. 

This recipe is described in grater detail below: 

First we calculate ylf '^^ and isf theo- 
retically using the full equations described in [20] . To ob- 
tain Ci we use CAMB [21, ^ with the WMAP 5-year best 
fit parameters and use the beam transfer functions from 
the WMAP team. We then obtain the connected piece 
using a modified version of the CMBFAST code [2^^. 
Plots of many of the quantities used for these calcula- 
tions can be found in Ref. [Hi. 
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with various assumed values 



of Tnl and ^nl to establish K-l ' ' and K-Y'^' . These are 
plotted in Figure [l] These curves will be compared with 
estimators derived from data to determine the magnitude 
of each statistic. Since we have two estimators, we can 
solve for the two unknowns tnl and jnl by fitting both 
estimators simultaneously. 

To calculate^ the estimators from data, used in the left- 
hand side of equations ([t]) and ([8| , we use both the raw 
and foreground-cleaned WMAP 5- Year Stokes I maps for 
V- and W-bands masked with the KQ75 mask ^. We use 
the Healpix library to analyze the maps. For this analysis 
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^ http ; //www. cfa.harvard.edu/ mzaldarr/CMBFAST/cmbfast.html 
^ see Smidt el al. 2009 for a similar calculation using the bispec- 

trum for more details. 1141 
* http : / /lambda. gsfc.nasa.gov/ 
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FIG. 2: The relation between the full estimators coming from 
data versus the Gaussian contributions. The green curve show 
the Gaussian contributions coming from averaging the estima- 
tors from the Gaussian maps. The red curve is the theoretical 
Gaussian piece calculated from Eq. ?? using the WMAP-5 
best-fit cosmology power spectrum. The error bars show two 
standard deviations from the Gaussian curves. These curves 
are from W band data. 



we only considered data out to ^max — 600. We correct 
for the KQ75 mask using a matrix Af;;/, based on the 
power spectrum of the mask, as described in [20| . 

Figure [l] shows the results for K?{^ and /Cf'^ for the V 
and W frequency bands extracted from the raw WMAP 
5- Year maps. In order to do proper statistics for our data 
fitting we create 250 simulated Gaussian maps of each 
frequency band with rigidc — 512. To obtain Gaussian 
maps we run the synfast routine of Healpix with an in-file 
representing the WMAP 5-year best-fit CMB anisotropy 
power spectrum and generate maps with information out 
to I — 600. We then use anafast, without employing 
an iteration scheme, masking with the KQlb mask, to 
produce aim's for the Gaussian maps out to I — 600. 
Obtaining estimators from these Gaussian maps allows us 
to uncover the Gaussian contribution to each estimator in 
addition to providing us information needed to calculate 
the error bars on our results. 

This whole process is computationally intensive. To 
calculate all theoretical estimators took nearly 8,000 
CPU hours. Furthermore, all the estimators from Gaus- 
sian and data maps combined took an additional 1600 
CPU hours. 

As previously discussed, the full trispectrum can be de- 
composed into both a Gaussian and non-Gaussian or con- 
nected piece. To make a measurement of non-Gaussianity 
we to subtract off the Gaussian piece from the full trispec- 
trum. Figure [2] shows the the relationship between the 
full trispectrum and the Gaussian piece. In this plot the 
Gaussian piece was calculated in two different ways as 
a sanity check. First, the Gaussian maps were averaged 
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FIG. 3: The 95% confidence levels for (?nl versus tnl- The red 
and orange represent the 68% and 95% intervals respectively 
for the combined V+W analysis. The light blue regions rep- 
resent the 95% confidence intervals for the V band analysis, 
and the light green regions are for the W band. 



Band 


W 


V 


V±W 


Raw 








ffNL 


4.7x10* ± 5.3x10^ 


4.6x10* ± 5.9x10'^ 


4.7x10* ± 3.9x10^ 


Tnl 


(1.63± 1.27)xl0'' 


(1.68± 1.31)xl0* 


(1.64±0.98)xl0* 


^NL 


7.4 ±7.3 


6.3 ±6.0 


11.1 ± 7.3 


FC 








5NL 


4.2x10* ± 5.3x10^ 


4.1x10* ±5.9x10^ 


4.2x10* ± 3.9x10^ 


Tnl 


(1.32± 1.27)xl0* 


(1.39± 1.31)xl0* 


(1.35±0.98)xl0* 


j4nl 


6.0 ±6.7 


5.2 ±5.7 


9.2 ±6.1 



TABLE I: Results for each frequency band to la. Values for 
Qnl, Tnl and j4nl on the top are for raw maps. The values 
on the bottom are for foreground clean maps. 



over. Second, the Gaussian piece of each estimator is 
calculated theoretically as described in|20j. 

After obtaining the theory, data and simulated curves 
we use the best fitting procedure described in [,14 where 
we minimize to fit tnl and gNL simultaneously. Our 
results are listed in Table |l| We see that (/nl and tnl 
are consistent with zero with 95% confidence level ranges 
-7.4 < ^nl/IO^ < 8.2 and -0.6 < tnl/IO"* < 3.3 for 
V±W-band in foreground-cleaned maps. The 95% confi- 
dence intervals of ^nl versus tnl are plotted in Figure [3] 
for each band. Furthermore, we find that the constraint 
-3 < Anl < 21.4. 

Conclusion. — This paper gives the first direct con- 
straints on Tnl and improves the previous indirect es- 



timate by four orders of magnitude. Our limit on tnl is 
close to a level where interesting constraints can started 
to be placed on models available in the literature for pri- 
mordial perturbations, such as due to cosmic strings[7]. 

Furthermore, this paper is the first to constrain Anl- 
This quantity is a model dependent relation between t^l 
and /nl , that for some models can be used to further con- 
strain the tensor-to-scaler-ratio. We uncover that that 
the data is consistent with Anl values between -3 and 21 
at 95% confidence. Moreover, models with Anl largely 
negitive, such as some arising from DBI infiation are dis- 
favored [7^. 
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